Strongly interacting two-dimensional Bose gases 
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We prepare and study strongly interacting two-dimensional Bose gases in the superfluid, the 
classical Berezinskii-Kosterlitz-Thouless (BKT) transition, and the vacuum-to-superfluid quantum 
critical regimes. A wide range of the two-body interaction strength 0.05 < g < 3 is covered by tuning 
the scattering length and by loading the sample into an optical lattice. Based on the equations of 
state measurements, we extract the coupling constants as well as critical thermodynamic quantities 
in different regimes. In the superfluid and the BKT transition regimes, the extracted coupling 
constants show significant downshifts from the mean-field and perturbation calculations when g 
approaches or exceeds 1. In the BKT and the quantum critical regimes, all measured thermodynamic 
quantities show logarithmic dependence on the interaction strength, a tendency confirmed by the 
extended classical-field and renormalization calculations. 
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Two-dimensional (2D) Bose gases are an intriguing 
system to study the interplay between quantum statis- 
tics, fluctuations and interaction. For non-interacting 
bosons in 2D, fluctuations prevail at finite temperatures 
and Bosc-Einstcin condensation occurs only at zero tem- 
perature. The presence of interaction can drastically 
change the picture. With repulsive interactions, fluctua- 
tions are reduced and superfluidity emerges at finite tem- 
perature via the Berezenskii-Kosterliz-Thouless (BKT) 
mechanism [l], In cold atoms, the BKT transition 
and the suppression of fluctuations are observed based 
on 2D gases in the weak interaction regimes 
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Extensive theoretical research on 2D Bose systems ad- 
dresses the role of interactions in the superfluid phase @- 
[l2j and near the BKT critical point [11 0. 

In the weak 

interaction regime, the classical 4 field theory [11 [l4| 
predicts the logarithmic corrections to the critical chem- 
ical potential /x c = k-QT{g/ir) ln(13.2/<?) and the criti- 
cal density n c = A^ B ln(380/g) for small two-body in- 
teraction strength g < 0.2. Here k-gT is the thermal 
energy and AdB is the thermal de Broglie wavelength. 
The classical-field predictions are consistent with weakly 
interacting 2D gas experiments Ell • 

In the superfluid ground state, Popov showed that 
the ground state chemical potential [i deviates from the 
mean-field result //mf = ti 2 gn/m logarithmically Q. 
Here m is the mass of the boson, n is the density, and 
2irh is the Planck constant. Defining the superfluid cou- 
pling constant as G — m/{h 2 n) where k = dn/d^i is the 
compressibility, we can summarize the perturbation ex- 
pansion result of G as Q 
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FIG. 1: (color online) Equations of state for 2D Bose 
gases and 2D lattice gases with 0.05 < g < 2.8. The 

filled circles represent measurements of 2D gases with g — 0.05 
(black), 0.15 (red), 0.24 (blue), 0.41 (green) and 0.66 (purple). 
The open circles represent measurements of 2D lattice gases 
with g = 0.45 (black), 0.85 (red), 1.2 (blue), 1.9 (green) and 
2.8 (purple). The blue shaded area is the superfluid regime, 
and the red shaded area corresponds to the BKT transition 
regime. The black dashed line fi = indicates where we 
evaluate the density and pressure for a vacuum-to-superfuid 
quantum critical gas. The inset compares the equations of 
state of a 2D gas and a 2D lattice gas with an almost identical 
g w 0.4. 



G = 



1 + Aging + Bg + Cg 2 \ng + Dg 2 



(1) 



where A = -I /4m [8], B = (In 4 - 27 - 2)/4tt M, C = 
— l/167r 2 0|, the value of D remains controversial |10l [TTIj . 
and 7 is Euler's constant. 

Beyond perturbation, calculations based on the renor- 
malized classical Ginzburg-Landau theory [ill [l7| at fi- 
nite temperature yield the result G = A recent 
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nonperturbative renormalization-group (NPRG) calcula- 
tion also provides complete thermodynamic calculations. 
Near the vacuum-to-superfluid quantum critical point, 
where the chemical potential /i = and the temperature 
T = 0, dimensionless pressure P is approximated to be 



P = 92 (e 



-( ff /9.1)W(9.X/ 9 ) 



), where 32(2;) = X^feli # /Ar is 



the Bose function. W(x) is the Lambert function satis- 
fying W{x)e w ^ = x, and the dimensionless density is 
n = - ln(l - e -(9/3-8W3.8/s)) 

In this letter, we extend our previous work on weakly- 
interacting 2D Bose gases Q into the strong interac- 
tion regime. We test the above theoretical predictions 
in different regimes, see Fig. 1, and our measurements 
show significant deviations from the mean-field theory 
as well as the logarithmic dependence on the interaction 
strength. 

A continuous evolution of a 2D quantum gas from the 
weak interaction {g <C 1) to the strong interaction {g > 1) 
regime is achieved by tuning the magnetic field near a 



Feshbach resonance 21], and by combining experiments 
with and without an optical lattice. Optical lattices en- 
hance the interaction strength by increasing the on-site 
density and the effective mass to*. The definition of g 
for a 2D gas (no lattice) is given in Refs. [22T - I241 ] . and 
for a 2D lattice gas given in Ref. (25|. Both definitions 
are mutually consistent, and can be connected to the 2D 
interaction strength g — 47r/| In na^L where a2D is the 
scattering length in two dimensions [26] . 

We start our experiment by preparing a degenerate 
Bose gas of cesium atoms in a two-dimensional optical 
trap [y, [25| . The radial and axial angular trap frequen- 
cies are (uj x ,u!y,u) z ) = 2tt x (8, 10, 1900) Hz. The sample 
contains 2x 10 4 atoms with temperature T = 13 ~ 20 nK, 
well below the excitation energy in the z— direction such 
that the sample is in the quasi-2D regime [22T |. We use 
a magnetic field to tune the atomic scattering length 
a — 40 ~ 580 ao <C l z near a Feshbach resonance, where 
ao is the Bohr radius and l z = 200 nm is the harmonic 
oscillator length in the z direction. The corresponding 
interaction strength is g = 0.05 ~ 0.77. 

To further enhance the interaction, we load the 2D gas 
into an optical lattice. A 2D square lattice is formed with 
a lattice constant of 532 nm, and the depth is set to be 
V = 7.1 -Er, = k-Q x 450 nK, where the tunneling energy is 
t = fcex 2.5 nK, the effective mass is to* = 2.9(1) to, 
is the recoil energy and fee is the Boltzmann constant. At 
this lattice depth, the system is far from the unity-filling 
Mott insulator phase and, for all interaction strengths 
we study, the ground state of the system remains in the 
superfluid phase. For 2D lattice gases, we can tune the 
interaction strength up to g — 2.8. 

To ensure thermal equilibrium, we prepare the gases at 
different interaction strengths by adiabatically ramping 
the magnetic field and the lattice potential. For all 2D 
gas experiments, we use a 200-ms magnetic field ramp 
which is slow compared to the time scale of the radial 
motion. For the 2D lattice experiments, we adopt an 
adiabatic lattice potential ramp of 400 ms [27| . The mag- 
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FIG. 2: (color online) Coupling constant Gsf of 
strongly interacting 2D superfluids. We determine Gsf 
by fitting the slope of the equations of state in the superfluid 
region for 2D gases (filled circles) and 2D lattice gases (open 
circles). Extensions of theoretical predictions into the strong 
interaction regime based on the mean- field theory (red line), 
the third-order perturbation expansion 7j| (green line), see 
Eq. (1), and the Ginzburg-Landau theory [17] (black line) are 
shown for comparison. The error bars are dominated by the 
uncertainty of the density calibration. 



netic field ramp is performed within the first 200 ms of 
the lattice ramp. For both the 2D gas and the 2D lattice 
gas, we monitor the subsequent density distribution for 
up to 200 ms after the ramp and observe no detectable 
dynamics and insignificant atom loss [2c| . 

We determine the equations of state by measuring in 
situ atomic density profiles based on absorption imaging 
with a high resolution objective (numerical aperture = 
0.5). Imaging aberrations are carefully characterized [29| . 
As a result, we achieve a spatial resolution of 1.0 /zm. The 
atomic density is calibrated by the number fluctuation 
of a normal gas Q. The measured density profiles are 
then converted into the equation of state n(/x, T, g) based 
on local density approximation [30j, where \i and T are 
determined by fitting the density tail @, H|| HH . Note 
that we define the zero of the chemical potential to be 
the energy of the lowest available single particle state in 
order to compare the equations of state of both 2D gases 
and 2D lattice gases. 

We plot the equations of state of 2D gases and 2D 
lattice gases in the dimensionless form n*(/z), where n* = 
nA*jg is the phase space density, X^ B = hj \J2-Kra*k^T is 
the thermal de Broglie wavelength and jl = fi/ksT is 
the dimensionless chemical potential }32j. For 2D gases, 
the effective mass is to* = to. Samples of the measured 
equations of state are shown in Fig. 1. In particular, two 
equations of state with a similar g s» 0.4, one from a 2D 
gas with a = 310 ao and one from a 2D lattice gas with 
a — 40 ao, are compared in the inset of Fig. 1. The overall 
matching behavior of the two equations of state justifies 
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FIG. 4: (color online) Pressure (a) and density (b) of a 
quantum critical gas at ft = 0. The measurements based 
on 2D gas (filled circles) and 2D lattice gas (open circles) at 
temperatures between T = 11 ~ 15 nK are compared with 
NPRG theory pl| (black line), mean-field theory for a 2D 
gas [35[ (blue line) and for a 2D lattice gas [36l ] (red line) at 
13 nK. 



FIG. 3: (color online) Critical coupling constant (a), 
scaled critical chemical potential (b) and scaled crit- 
ical density (c) determined in the BKT transition 
regime. By overlapping all scaled equations of state in the 
transition regime, shown in the panel (a) inset, critical pa- 
rameters are determined from Eq. J2]). The results from 2D 
gases (filled circles) and 2D lattice gases (open circles) are 
compared to the predictions from the mean-field theory (red 
line), the perturbation theory [?J (green line), the classical- 
field theory [l3|] (blue line) and the Ginzburg-Landau theory 
[13 (black line). 

our use of optical lattices to enhance the interaction. The 
small discrepancy near ft « will be discussed below. 

In the superfluid regime, we extract the coupling con- 
stant Gsf = m* / (H 2 k) by evaluating the superfluid com- 
pressibility k — dn/dfi, see Fig. 2. The coupling con- 
stants show significant downshifts from the mean-field 
prediction when the system enters the strong interaction 
regime. Similar tendency is also shown in the Ginzburg- 
Landau calculation [17j . as well as in a recent work [33j ] 
which includes effective three-body interactions. 

In the BKT transition regime, we use the universal 
critical behavior of the equations of state to determine 
the critical parameters 0, l2fjj . By rescaling and over- 
lapping [34j ] all the equations of state in the transition 
regime according to 



h ._ K = H{ !L_?i )i (2) 

we obtain the critical phase space density h* , the critical 
chemical potential ft c and the critical coupling constant 
G c , see Fig. 3. Remarkably, equations of state of all 2D 
gas and 2D lattice gas measurements overlap excellently, 
see Fig. 3 (a) inset. 

The extracted critical coupling constants G c are con- 
sistently lower than the mean- field values G — g. On 
the other hand, the extracted scaled critical chemical 
potentials ft c and scaled critical densities h* follow the 
logarithmic dependence on g predicted by the classical- 
field calculations [U, HH E3]- Our results confirm the 
crucial role of interactions in 2D Bose gases and suggest 
that the extensions of the above theories into the strong 
interaction regime capture the general behavior of the 
thermodynamic quantities. 

Finally, we investigate the pressure and density in the 
quantum critical regime. In the lattice, atoms reach 
the vacuum-to-superfuid quantum critical regime when 
the chemical potential matches the lowest single particle 
state, and when the thermal energy is below the ground 
band bandwidth (25|. We extend the definition of quan- 
tum criticality to 2D gases [23]. To determine the pres- 
sure, we integrate the density over the chemical poten- 
tial, Pq = f_ n*dft. The extracted P * in 2D gases and 
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2D lattice gases are compared with the mean-field and 
NPRG calculations [T^], see Fig. 4 (a). Here, we observe 
overall agreement between experiment and theories. For 
lattice gases, in particular, the slightly higher even in 
the weak interaction regime is discussed in Ref. [19| as the 
result of finite temperature effect. The densities in the 
quantum critical regime Uq also show the expected log- 
arithmic dependence on the interaction strength. Here, 
we observe small systematic deviations from the theories. 

To conclude, we report the preparation and thermody- 
namic study of stable strongly interacting 2D gases. Di- 
mensionless coupling constant g as high as 2.8 is reached 
by Feshbach tuning in an optical lattice. In the strong 
interaction regime, coupling constants show clear devi- 
ations from the mean-field theory. Other thermody- 
namic quantities in the classical and quantum critical 



regimes show strong dependence on g and can be cap- 
tured well by extensions of the classical-field theories and 
the NPRG calculation. Our results provide new insight 
into the crucial role of interactions in the thermodynam- 
ics of 2D gases. Further enhancement of the interaction 
strength can potentially lead to crystallization of the 2D 
gas [33] . Investigation on the fluctuation and correlation 
of strongly interacting 2D gases will be reported else- 
where. 
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